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ABSTRACT 

Meteoritic chondrules were formed in the early Solar System by brief heating of silicate dust to melting tem- 
peratures. Some highly refractory grains (Type B calcium-aluminum rich inclusions, CAIs) also show signs of 
transient heating. A similar process may occur in other protoplanetary disks, as evidenced by observations of 
spectra characteristic of crystalline silicates. One possible environment for this process is the turbulent magne- 
tohydrodynamic flow thought to drive accretion in these disks. Such flows generally form thin current sheets, 
which are sites of magnetic reconnection and dissipate the magnetic fields amplified by a disk dynamo. We 
suggest that it is possible to heat precursor grains for chondrules and other high-temperature minerals in current 
sheets that have been concentrated by our recently described short-circuit instability. We extend our work on 
this process by including the effects of radiative cooling, taking into account the temperature dependence of the 
opacity; and by examining current sheet geometry in three-dimensional, global models of magnetorotational in- 
stability. We find that temperatures above 1600 K can be reached for favorable parameters that match the ideal 
global models. This mechanism could provide an efficient means of tapping the gravitational potential energy 
of the protoplanetary disk to heat grains strongly enough to form high-temperature minerals. The volume-filling 
nature of turbulent magnetic reconnection is compatible with constraints from chondrule-matrix complemen- 
tarity, chondrule-chondrule complementarity, the occurrence of igneous rims, and compound chondrules. The 
same short-circuit mechanism may perform other high-temperature mineral processing in protoplanetary disks 
such as the production of crystalline silicates and CAIs. 

Subject headings: instabilities — magnetic reconnection — magnetohydrodynamics — plasmas — protoplan- 
etary disks 



1. INTRODUCTION 

The meteoritical record indicates that chemical and min- 
eralogical evolution of protoplanetary disks depends on lo- 
cal heating processes. The ubiquity of chondrules in mete- 
orites, along with the relation between chondrules and matrix 
(including chondrule-matrix elemental complementarity), ap- 
pears to demand a heating mechanism that acts locally, near 
the assembly point of chondritic meteorites. Type B calcium- 
aluminum rich inclusions (CAIs), found in CV chondrites, 
also appear to have experienced transient heating, though un- 
der different conditions. 

The basic challenge in producing chondrules is heating 
precursor grains with radii ~ 1 mm from below 1000 K 
to ~ 1800 K, and then cooling them at rates of 100- 
1000 K hour~^, significantly sl ower than expected for free 
radiation to a cold backgrou nd ([Hewins & Radomsky 19901 
ILofgren & Lanier 1990; Connolly et al. 1998; , Scott & Krot| 
|2005 ; Ebel 2006 ). It further appears that heating to - 2000 K, 
followed by cooHng to 1500 - 1800 K at ^ 5000 K hour" -^ 
is necessary to form barred textures ( [Connolly et al.|[20Q6l ). 
Particularly low starting temperatures around 650 K may be 
required to retain the sodium and other volatile elements that 



cmcnally @ nbi.dk, ahubbard @ amnh.org^ mordecai @ amnh. org,] 

r debel@amnh.org, p.dalessio@crya.unam.nix 

^ Niels Bohr International Academy, Niels Bohr Institute, Copenhagen, 
Denmark 

^ Department of Astrophysics, American Museum of Natural History, 
New York, NY 10024-5192, USA 

^ Department of Astronomy, Columbia University, New York, NY 
10027, USA 

^ Department of Earth and Planetary Sciences, American Museum of 
Natural History, New York, NY 10024-5192, USA 

^ Centro de Radioastronomia y Astrofisica, Universidad Nacional 
Autonoma de Mexico, 58089 MoreHa, MICH, Mexico 



are found in some chondrules ( [Alexander et al.|2008] ), as these 
volatile elements will be released when the chondrules are 
melted. Regardless of the starting temperature, strong clus- 
tering of dust grains and rapid processing is required to avoid 
depletion ( Alexander et a l. 2008 ). 

Examination of the chondrules and matrix in several chon- 
dritic meteorites has shown that the two constituents appear 
complementary in that, though the elemental makeup of each 
varies substantially, the bulk chemical compositions of chon- 
drites made up o f both always reflects the eleme ntal ratios of 
the sun (Korna cki & Wood|[r984l [Bland et al.| 2005 ; HezeT[ 



& 



Palinej|2010|. Indeed, chondrules with high variability 
may themselves be complementary to each other. 



may tnemseives be complementary 
combining into parent asteroids with solar Fe/Si ratios (Ebel] 
et al. 2 008 ). This suggests that chondrules must have formed 
from local heating of matrix material. Grains from other 
stars (presolar grains) and insoluble organic matter found 
in the matrices of most primitive chondrites co uld no t have 
survived the heating experienced by chondrules ( [Mendybaev[ 
[eFal. 2002 i Connolly et al. 2006 ), establishing that the heating 
events were intermittent and did not process all available ma- 
terial. Many chondrules experienced multiple heating events, 
as shown by the occurrence of igneous rims ( Jones 2012 ), and 
collisions with partiall y molten and solid droplets to produce 
compound chondrules ( [Ciesla et al.|2004[ ). This supports the 
hypothesis that the chondritic meteorites are formed soon af- 
ter chondrule melting, out of a dense cluster of dust grains 
only a fraction of which are melted. Such a scenario requires 
a highly localized heating source. 

Type B CAIs are remelted highly refractory inclusions. 
They show evaporation of Mg and Si due to melting and 



evaporation in reheating ev ents af ter condensation ([Grossman 
et al.|2000t|Richter et al.|2002HShalw"& Young|2007[ ). Re 
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cent evidence sugg ests that the epochs of CAI and chondrule 
formation overla p ( [Moynier et al.|20Q7[[Yin et al.|20Q9[|Con 
nelly et al |20T? ). 



A possible energy source for this localized heating is the 
magnetorotational instability (MRI), which taps the orbital 
energy in the disk's differential rotation to drive magnetized 
turbulence, transporting ang ular momentum outward and al- 



lowing accretion to occur (|Velikhov||1959t 

)9irThe MRI cf] 



jChandrasekhar 
riven turbulence 



T96Ql|Balbus & Hawley| [19911 ) The 
dissipates energy as heat. However, the heating occurs in- 
termittently, not uniformly. MHD turbulence quite gener- 



ally f orms current sheets (Parker 1972, 1994; jCowley et al. 
1997| ), that dissipate energy faster than the average rate. Mag 



netic reconnection can occur in such current sheets, enhancing 
their local heating rate in small regions throughout the cur- 
rent sheets (For a review, see Yamada et al. 2010). Hub bard] 
|et al.| ('2012 , hereafter Paper I) showed that if the ionization 
depends strongly on the temperature, the dropping resistivity 
within the region heated by a current sheet can drive an insta- 
bility effectively resulting in a short circuit, leading to strong 
localized heating. 

One of the main diffusive processes limiting where MRI 
occurs in a protoplanetary disk is Ohmic resistivity, which 
is in turn a function of the degree of ionization of the gas. 
Over much of the disk, estimates suggest that the ionization 
near the midplane is too low for the flow to be unstable to 
the MRI. It is hence expected that magnetically inactive dead 
zones exist in these den se, cool, and almost entirely neutral 
regions ( Gammie'1996). The secular movement of these dead 
zones through the disk may temporarily suppress the MRI, 
and thus transient heating, potentially explaining the broad 
range of observed chondrule ages. 

The most common treatment of ionization of a protoplan- 
etary disk includes two main components, thermal ionization 
of alkali metals (most importantly potassium), and nonther- 
mal ionization from radiation (stellar u ltraviolet light an d x- 
rays, cosmic rays, and radionuclides; Armitage |2011| ). It 
has also been proposed that non-ideal MHD effects beyond 
Ohmic res istivity limit MRI in various regions of the disk 
(|Blaes & Balbus 1994; Bai 2011a), and that the smallest dust 
grains can themselves be the dominant charge carriers in the 
outer disk (Bai2011b|. 

This letter first presents an analysis of typical current sheet 
geometries in MRI driven turbulence, and then presents one- 
dimensional models including radiative transfer, that demon- 
strate that the short-circuit instability first described in Paper I 
can produce sufficiently high temperatures in these current 
sheets to implicate them in the formation of high-temperature 
minerals. This instability is controlled by the rapid drop in re- 
sistivity with increasing temperature that results from thermal 
ionization of potassium at T > 1000 K. 

Other mechanisms for melting chondrule precursor s by 
means of magnetic reconnection have been proposed. Sonett 
([1979 ) suggested that beams of relativistic electrons emerg- 
ing from reconnection regions above or below the dense disk 
could heat chondrule precursors. The reconnecting magnetic 
fields were hypothesized to originate either from the Sun, or 
a disk dynamo. Along with a model for the inte rmittency 
of reconnection in a protoplanetary disk. King & Pringle 
(12010) proposed that magnetic reconnection regions could 
drive shocks strong enough to melt chondrule precursors. For- 
mation of chondrules by am bipolar diffus ion dominated re- 
connection was explored by Joung etaL] (|2004j). However, 
that mechanism is strongly suppressed by even a small guide 
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Figure 1. A typical reconnection region found in the global MRI model. The 
gray isosurface outlines high current density. Upper: Streamlines show mag- 
netic field, colored by azimuthal component. The current sheet is bordered 
by oppositely directed azimuthal flux tubes. Lower Left: The same viewpoint 
as the upper panel, but the streamlines are now colored by plasma /3. The 
value of ^ in the flux tubes bordering the current sheet is ~ 1, and at least 
~ 500 in the current sheet. Lower Right: Same quantities as the upper panel, 
with the viewpoint rotated to view down the short axis of the current sheet. 
The directions of the magnetic field in the bounding flux tubes can be seen 
to be nearly completely opposed. Field lines originating from points inside 
the current sheet where the field is weak wander significantly away from the 
directions of the bounding flux tubes. 



dMaron et al |20 12] [McNally et al.|2QT2| ), to perform a 
al, unstratified simulation of the MRI in a protoplanetary 



field in the reconnection region ( Heitsch & Zwei bel 2003 )^ 

The electrica l discharge mechanis m of [Inuts uka & S ano| 
( |2005| ) and Mur anushi et al.| ( 2012[ ) may mteract with the 
short-circuit instability, being triggered when the current sheet 
is sufficiently intensified. 

In Section |2] we describe the bounding geometries of cur- 
rent sheets analyzed in a global zoom-in simulation of MRI 
driven turbulence. In Section|3]we demonstrate that the short- 
circuit instability can heat these current sheets to chondrule- 
forming temperatures, and that the temperature variation of 
opacity can lead to a self-regulated form of the instability. We 
discuss the results, implications, and what further work must 
be done to establish short circuits as a heating mechanism in 
protoplanetary disks in Section [4] 

2. CURRENT SHEET GEOMETRIES 

We used Phurbas, an adaptive, Lagrangian, m eshless, MHD 
code \ ~ ' 
globaf 

disk, in order to examine current sheet strength and geome- 
try. An MRI-unstable cylindrical Keplerian shear flow with 
maximum radius 2.6 AU is placed in the center of a vertically 
periodic computational volume 40 x 40 x 0.4 AU. An isother- 
mal equation of state is used, with sound speed corresponding 
to a Keplerian disk with aspect ratio of 0.1, typical of proto- 
planetary disks. Initially, a vertical magnetic field is imposed, 
with a strength such that the critical MRI wavelength is one 
quarter of the vertical height of the volume. 

Once the MRI-driven flow is reasonably saturated (as indi- 
cated by the evolution of the magnetic tilt angle), two nested, 
crescent- shaped, more highly refined regions that rotate with 
the flow are added, centered at a radius 1 AU, and extending 
~ 7r/2 radians. These achieve resolution 1.7 and 3.4 times 
greater than the base flow. The Phurbas resolution parameter 
A = 10~^ AU in the lowest resolution region, and decreases 
to 2.8 X 10~^ AU in the highest resolution region. This is the 
region where current sheets were examined . Further details 
are provided in Chapter 6 of McNally '('20 12 ). 

In Figure [T] we show a typical high current-density region. 
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Figure 2. Heating of a current sheet bounded by (3 = 1.5, starting at To = 850 K. The current sheet narrows though the short-circuit instability as the resistivity 
7] drops. The peak temperature reached is T ~ 1650 K. Shown are magnetic field B, current density J, temperature T, resistivity 77, potassium ionization 
fraction Xe/a, inverse mean photon free path density p, velocity Vx, and total pressure Ptot, with respect to spatial position x and time t. The results are 
plotted mirrored about x = following the boundary conditions, and only the central section of the spatial domain is shown. 

which takes the form of an approximately two-dimensional 
current sheet sandwiched by large, nearly azimuthal flux 
tubes, whose magnetic field approaches the maximum field 
at their orbital position. The bottom panels show that there is 
very little perpendicular guide field passing through the cen- 
ter of the sheet, as the magnetic fields on the opposite sides 
are nearly perfectly opposed. In subsonic MHD turbulence, 
the expectation is that the peak currents occur at the tail of a 
turbulent cascade: spatially intermittent, short in duration and 
length scale, and randomly oriented. Peak magnetic fields 
generated by the MRI are, however, near equipartition, which 
may explain the large scale high current density regions. 

We find that the current sheets retain the same general struc- 
ture as the resolution is increased, while becoming thinner and 
more concentrated. The total, gas plus magnetic, pressure is 
approximately constant across these sheets, so in this locally 
isothermal simulation the sheets are regions of high density. 
From these results we conclude that a Harris-like current sheet 



geometry ( [Harris [1962 ) as invoked in Paper I, bounded by 
magnetic fields with magnetic pressure comparable to thermal 
pressure, is well motivated for the study of Ohmic heating in 
protoplanetary disks. This configuration uses a tanh profile 
for the magnetic field across the current sheet, and the total 
pressure is set constant by satisfying an adiabatic-hydrostatic 
condition with the temperature and density. 

3. ONE DIMENSIONAL MODELS 

We extended the implicit, one-dimensional, MHD code 
used in Paper I, which uses sixth-order finite differences on a 



logarithmic grid, and the CVODE package ( Hindmarsh et al. 
12005!) for time integration. The extensions include radiative 
transfer, and a full Saha-equation based treatment of the ion- 
ization of potassium (with abundance relative to hydrogen of 
a = 10 and background nonthermal ionization. The non- 
thermal ionization is parameterized following equation (14) of 
Froman g et al. ( 2002 ), with the ionization rate ( = 10~^^ s~^, 
an order of magnitude higher than the canonical value from 
cosmic ray ionization ( = 10~^^ s~^. This nonthermal ion- 
ization component is used to maintain low enough resistivity 
to hold our artificial current sheets in place, in the absence of 
an active turbulent flow. This slightly increases the onset tem- 
perature of the instability and slightly weakens its growth rate, 
so that our model establishes a lower limit to susceptibility to 
the instability. 

The code used in this Letter includes a radiative heating 
and cooling term calculated with the ray tracing scheme pre- 
sented in Heinemann et al.| (|2006 ). This radiative cooling is 
along one dimension, as the global simulations have shown 
that the transverse dimensions of the current sheets are much 
shorter than the parallel dimensions. This allows local Ohmic 
heating at the center of a current sheet to heat its surround- 
ings. With this additional physics, we have found it expedi- 
ent to add a grid- scaled sixth-order hyperdiffusion to the mo- 
mentum and temperature fields. Hyperdiffusions are scaled 
locally to the grid spacing. A logarithmic grid is specified 
by mapping a linear grid g which ranges from to Lg to the 
physical coordinates x by x{g) = fgLg{ex.-p{g/{fgLg)) — 1) 
where fg is a factor which controls the grid stretching. The 
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Figure 3. A thinner current sheet with width 10^° cm, with all other parameters identical to FigurejS] in which the potassium ionization saturates. The results 
are plotted mirrored about x = 0, and only the central section of the spatial domain is shown. 



diffusion operator for a variable / is then specified directly 
as df/dt = DgAg'^d^ f /dg^, where is the grid spacing 
increment in g, so that the diffusion scales with the grid. We 
also artificially smooth the variation in opacity through the ad- 
dition of a finite relaxation time physically motivated by the 
finite vaporization time of the dust grains. Rosseland mean 
opacity is evolved hy dhz/dt = {l/t^){hZr — i^) where hz is 
the opacity used in radiative transfer, hh'P IS the value of opac- 
ity derived from the opacity tables for the gas-dust mixture, 
and tf^ is a time constant. These alterations maintain the in- 
tegrability of the system of equations, and we have checked 
that the numerical smoothing coefficients are small enough to 
not significantly affect the results presented here. Our initial 
configuration is the same as in Paper I. 

The Rosseland mea n opacities Kr t hat we use are an up- 
dated version of the |D'Alessio et al.| (|2Q01) model. Impor- 
tantly, they incorporate a distribution oT dust particles com- 
prised of ice, organic, and silicate components. The sublima- 
tion temperatures of the different components are taken to be 
independent of the grain size, so that the opacity has a sharp 
transition downwards when these temperatures are reached. 
This model includes: graphite and "astronomical" silicates, 
with optical constants a nd fractional abundances from Wein- 
, and water ice with optical properties 



partner & Draine|(|20Ql 
from I Warrenfp98?| r^ 



The grain size distribution is n(r) ex 
^-3.5 between Tmin = 0.005 /im and Vmax = 1 mm, appro- 
priate for chondrule forming regions. The opacity is calcu- 
lated using Mie theory and assuming grains are spheres. The 
sublimatio n temperatures for sil icates and water ice grains are 
taken from|Pollack et al.|(|1994|), and for graphite, we assume 



Tsub = 1200 K. For the gas component, which is important 
when dust sublimates, we assume LTE. The dust-to-gas mass 
ratio for silicates is 4 x 10~^, for organics, 2.5 x 10~^, and 
for water ice, 4.7 x 10~^, so in total it is 1.21 x 10~^. Our 
model assumes a constant mean molecular mass /i = 2.32, 
which is not strictly consistent with the vaporization of solids. 
However, we stop calculations before the dissociation of H2 
introduces a significant effect. 

As discussed in Paper I, the short-circuit instability requires 
a strong negative dependence of the resistivity on tempera- 
ture. We start our simulations with T = Tq = 850 K, where 
the ionization already has a rapid dependence on temperature 
even outside the current sheet. 

In Figure |2] we show the evolution of a current sheet 
model with an initial ratio of thermal to magnetic pressure 
at the box edge of (3 = 1.5, background temperature Tq = 
850 K, and background density po = 10~^ g cm~^ with 
B oc tanh(x/5 x 10^^ cm). In this configuration, the short- 
est wavelength for MRI, based on the length scale where 
the magnetic Reynolds number is unity, is on the order of 
1.7x10^^ cm. The grid had 1000 points, = 5x10^ cm and 
fg = 0.09. Values of Dg = 10^ cm^ in the temperature 
equation, and Dg = 10^^ cm^ in the momentum equation 
were used. The opacity was smoothed with = 2 x 10^ s. 
Notice in particular the small width of the narrowed current 
sheet. As the opacity drops, the inverse mean free path of 
photons X (plotted in units of cm~^ in Figure [2]) decreases 
dramatically, and an optically thin bubble is formed around 
the current sheet as the temperature where silicates become 
liquid is passed. The low opacity bubble has small internal 
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temperature gradients, and hence smaller resistivity gradients. 
As the resistivity gradient weakens, the instability shuts off. 
In this way the instability self-regulates. For the parameters 
in Figurej2| it produces heated regions with temperatures that 
only slightly exceed the melting temperature of the silicate 
grains that dominate the opacity, the largest of which produce 
chondrules; at this point, the instability is saturated by the 
increased radiative cooling. This is the most obvious configu- 
ration leading to chondrule formation. 

If the current sheet is initially thinner, the heating rate be- 
comes fast enough that the heating is not affected by the 
drop in opacity. Figure |3] shows the same configuration as 
before, except with a factor of five thinner initial current 
sheet B ex tanh(x/10^^ cm). The grid had 800 points, 
Lg = 10^ cm and fg = 0.09. Values of Dg = 10^ cm^ s"^ in 
the temperature equation, and Dg = 10^^ cm^ in the mo- 
mentum equation were used. The opacity was smoothed with 
t/^ = 1 X 10^ s. In this case, the short circuit instability nar- 
rows the current sheet until the potassium ionization fraction 
(Xe/a) reaches unity and the resistivity gradient disappears. 
The current sheet continues to heat as it dissipates, although 
we stop our calculation as it breaks down because we do not 
include the effects of H2 dissociation and other changes to the 
effective equation of state above 2000 K. 

4. DISCUSSION AND CONCLUSIONS 

The models in this paper are far from a full treatment. How- 
ever, they demonstrate that the short-circuit instability is a 
plausible candidate for chondrule formation as shown by the 
model in Figure [2j which starts within the range of bound- 
ary conditions for current sheets seen in our global models of 
protoplanetary disks, and climbs to chondrule melting tem- 
peratures within hundreds of hours. For thinner initial current 
sheets, the instability can reach temperatures high enough to 
fully ionize potassium, which is also high enough to produce 
barred olivine textures. This behavior is also compatible with 
the remelting of Type B CAIs. As the formation of the first 
chondrules and CAIs overlaps in tim e ( Moynier et al.||2007 



Yin et al.|2009l [Connelly et al.|2012| ), it is plausible that the 



same mechanism is responsible for processing Type B CAIs 
and the oldest chondrules. 

We have also noted that the drop in opacity associated with 
the destruction of organic components, at 1200 K in our opac- 
ity model, can produce a behavior similar to that in Figure [2] 
but at temperatures below the melting point of silicates. The 
temperatures produced in this type of event would be appro- 
priate for quickly annealing amorphous silicate dust to pro- 
duce crystalline silicates. 

Chondrule formation by the melting of precursor grains in 
reconnection regions subject to the short-circuit instability of 
Paper I seems compatible with constraints from chondrule- 
matrix complementarity, chondrule-chondrule complementar- 
ity, and multiple heating listed in Section[T} as it is a local pro- 
cess that can occur many times to the same body of meteorite 
parent body precursor material. 

As the MRI generates magnetic fields near equipartition 
with the gas thermal energy, compressive heating is a sub- 
stantial contributor to the temperature. Our initial condition, 
as in Paper I, is chosen to be adiabatic-hydrostatic, which is 
the gentlest configuration we have tested. Determining the 
strength of the compressive heating, and its ability to meet the 
threshold temperature for the short-circuit instability, will re- 
quire models that self-consistently form the current sheets in a 
large scale simulation that includes the action of temperature- 



dependent resistivity. Such simulations will also allow de- 
termination of the correct physical length scale for the current 
sheet thickness at the onset of the instability. This length scale 
will strongly affect the heating and cooling times. 

Furthermore, multidimensional studies of the current sheets 
are needed, to establish how they might break up (due to plas- 
moid instabilities or buckling of the current sheet) and hence 
the size and lifetime of the heated regions, which will be im- 
portant for determining the cooling rates of chondrules in this 
scenario. Finally, we note that the temperatures explored here 
have been discussed in reference only to chondrule and CAI 
formation. However, given more complete ionization models, 
other high temperature mineral processing, such as the crys- 
tallization of silicates, destruction of organics, or sublimation 
of ices, may occur through the same mechanism. 
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